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EXECUTIVE  SUMMARY 


This  report  describes  dynamic  combat  models  that  reflect  the  effect  of 

information  flows  together  with  attrition  capability  upon  combat  progress  and 

outcome.  Command  and  Control  assets  for  each  participant  are  modeled  as  endowed 
.  2 

with  the  capacity  to  guide  combat;  C|3is  also  vulnerable  in  that  it  may  be  delib¬ 
erately  targetted  and  reduced  In  effectiveness.  Physical  attrition  is  modeled 
first  by  a  deterministic  rate  process  (Lanchesterian  in  nature),  secondly  by  a 
stochastic  process  related  to  the  first  by  ideas  related  to  those  of  stochastic 
difference  and  differential  equations. 

The  models  are  best  exercised  and  explored  on  an  interactive  computer  display 
A  FORTRAN  program  exists  for  this  purpose,  with  displays  of  hypothetical  "histori¬ 
cal*  combat  outcomes  now  appearing  in  tabular  form.  Graphical  displays  will  be 
provided  in  future  work. 

It  seems  likely  that  this  model  simulation  can  be  developed  into  a  gaming 
tool,  very  conveniently  playable  by  two  persons  who  can  elect  various  strategies 
for  force  allocation,  play  the  game,  and  learn  from  the  results.  An  elaboration 
of  the  model  may  serve  as  a  means  for  assessing  the  importance  of  increased 
effectiveness  of  equipment,  either  with  respect  to  firing  rate  and  accuracy  or 
time  for  information  flow.  , 
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MODELS  OF  CONFLICT,  WITH  EXPLICIT  REPRESENTATION  OF  COMMAND  AND  CONTROL 
CAPABILITIES  AND  VULNERABILITIES 

Donald  P.  Gaver 


1.  INTRODUCTION 

Many,  if  not  most,  conflicts  between  opposing  forces  R  and  B  are  conducted 

3 

under  some  form  of  Command,  Control  and  Communications  (C  )  establishment  super¬ 
vision.  Yet  few,  if  any,  simple  analytical  models  seem  to  attempt  portrayal  of 

3 

the  relationship  between  C  and  combat  effectiveness.  This  report  suggests 

3 

models  involving  C  capabilities  and  vulnerabilities,  and  indicates  the  manner 
in  which  the  models  suggested  may  be  utilized  in  a  gaming  context. 

Models  constructed  in  the  present  manner  were  proposed  by  Gaver  and  Tonguij 
(1979).  In  that  study  the  opposing  forces  were  each  split  into  two  groups  char¬ 
acterized  by  their  respective  information  states:  those  "in  the  know,"  or 
capable  of  engaging  in  efficient  attrition  activities,  and  those  "in  the  dark," 
and  capable  only  of  less  appropriate  action,  or  none  at  all.  It  was  illustrated 
by  Tongu$  (1980)  that  a  capability  for  quick  transition  from  one  information 
state  to  another  could  sometimes  outbalance  raw  physical  capability,  such  as 
firing  rate  and  single-shot  kill  probability,  thus  acting  as  a  "force  multiplier. 

Once  the  potential  of  the  C  component  is  recognized,  the  latter  also 
becomes  a  potential  target.  It  is,  therefore,  of  interest  to  incorporate  the 
C  component  explicitly  into  attrition-type  models,  and  then  to  exercise  the 
models  so  as  to  expose  vulnerabilities  and  possibly  suggest  sensible  doctrine. 

We  make  a  stab  at  this  program  here  in  a  highly  simplified  manner,  feeling  that 
informative  elaborations  may  well  be  suggested  after  an  initial  look.  The  idea 
is  to  stray  in  a  gingerly  fashion  into  the  area  between  classical  Lanchester 
combat  theory  and  the  extensive  and  elaborate  terrain  of  the  modern  wargame. 


2.  MODELS  OF  POSITION  DEFENSE 

2.1.  Modal  I:  Defense  of  a  Stronghold  (or  Bastion  or  Beachhead) 

An  initial  model  for  this  situation  was  presented  by  Gaver  and  Tongue  (1979); 
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see  Tongue  (1979),  henceforth  call  this  the  GT  Model.  Suppose  an  R-force  of  size 

9 

R  attacks  a  bastion  (e.g.  beachhead  or  defended  position)  held  by  B.  The  sig¬ 
nificance  is  that  B  is  advantageously  located,  and  in  GT  was  assumed  to  suffer 
no  casualties  initially,  while  R  is  exposed  and  vulnerable  and  can  only  succeed 
by  (a)  surprise,  or  (b)  B’s  inability  to  critically  diminish  R  before  being 
overrun.  Suppose  that  R's  speed  of  advance  is  (nearly)  constant,  and  that  R 
wins  if  R(tg)  >_  kB,  where  RUq)  is  R's  force  size  when  the  stronghold  is 
reached  at  time  tQ;  k  represents  R's  necessary  advantage  over  B  at  final 
stages  or  "hand-to-hand"  in  order  to  win. 

In  the  GT  model,  (and  here)  it  was  (and  is)  assumed  that  B's  divide  into 
two  combat  groups:  one  in  number  Bu ( t ) ,  e.g.  the  number  of  those  able  to  fire 
in  an  ineffective,  specifically  unaimed,  manner,  and  another  of  size  B  (t), 

a 

e.g.  the  number  of  those  capable  of  firing  in  a  more  effective,  specifically 
aimed  way.  All  sorts  of  refinements  are  possible,  but  for  the  present  two 
information  states  are  sufficient. 

In  GT  it  was  assumed  that  dynamic  transition  between  B, (t)  and  B  (t) 

li  a 

3 

occurred:  the  rate  of  transition  B  (t)  -*•  B  (t)  measured  the  power  of  the  C 

U  a 

system.  However,  no  attempt  was  made  to  represent  that  system  as  an  explicit 
entity,  itself  being  vulnerable  and  hence  an  inviting  target.  In  the  present 

3 

model  assume  that  the  C  assets  of  B  are  vulnerable  to  R,  and  that  the  rate 

3 

of  transition  from  B  (t)  to  B,(t)  is  made  possible  by  the  C  force  B  (t);  the 

U  a  C 

latter's  effectiveness  can  in  turn  be  effected  by  R's  actions.  In  other  words,  R 
can  attempt  to,  or  inadvertently,  target  Bc(t)  --  and  only  Bc ( t )  in  the  present 
model  --  using  the  force  RCCU)  assigned  for  that  purpose.  To  the  extent  that 


Rcc(t),  called  the  Counter-C3  (C-C3)  Force,  is  effective,  the  C3  capability  of  B, 

namely  Bc(t),  is  reduced  by  temporary  suppression  or  outright  destruction.  Such 

reduction  in  turn  adversely  affects  the  quality  of  B's  response  to  R's  attack. 

Since  in  this  model  the  effort  (in  terms  of  force  size)  allocated  to  C-C  activity 

is  removed  from  the  Red  active  list,  thus  diminishing  R^  and  hence  the  number 

available  to  encounter  B  once  the  bastion  is  reached,  there  is  clearly  a  trade- 

3 

off  opportunity  for  R.  Too  large  an  active  force  R^  at  the  expense  of  C-C  , 

R  ,  allows  extensive  attrition  by  B  on  R.  free  of  charge:  B  (t)  may  work  at 

tv  M  v 
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full  effectiveness.  On  the  other  hand,  too  large  a  C-C  force  obviously  penal¬ 
izes  the  attacking  force,  RA  .  Similar  choices  exist  for  B.  Only  by 
setting  down  a  quantitative  representation  of  the  combat  dynamics  and  studying 
its  implications  numerically  can  informed  intuition  be  developed  that  may  lead 
to  a  wise  trade-off. 

Lanchester-Style  Equations 

Here  are  some  SDecific  Lanchester-style  differential  equations  that 
represent  the  dynamics  described  above  in  words. 

3  3 

•  Blue's  C  Interactions;  Red's  Counter-C  . 


u'  '  _ 
~3t~  - 

-  CM(Bc(t).  Bu(t» 

(2.1) 

dB.(t) 

a 

dt 

cua<Bc(t>'  Bu<*» 

(2.2) 

dB.(t) 

dt 

Ccc(Bc(t>-  Rcc(t>> 

(2.3) 

These  equations  represent  the  rate  of  change  of  the  B-force  segments  with  time. 
They  describe  only  C3-related  activities,  since  physical  change  (attrition)  of 
Bu  and  Ba  has  not  been  allowed  in  the  present  scenario.  Expression  (2.1)  states 
that  the  rate  of  change  reduction  of  the  Blue  ineffective  (unaimed)  force  depends 


3 


3 

through  the  function  Cua(.  ,  .)  upon  the  capacity  of  the  Blue  C  activity,  measured 

by  Bc(t),  and  upon  the  number  of  ineffective  forces,  Bu(t) ,  awaiting  conversion 

to  the  effective  state.  Notice  that  this  expression  is  left  general;  Cua(.  ,  .) 

3 

can  be  specified  at  will,  and  must  represent  the  general  features  of  the  C 
activity,  including  sensor  performance  and  output  analysis  as  well  as  communication. 
One  simple,  tentative,  but  specific  representation  might  be 

Cua<Bc(t)-  Bu(t»  -  6uaBc(t)  "„<*>• 

e  being  a  positive  constant.  This  expresses  the  appealing  intuition  that  the 

lid 
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rate  of  transfer  of  B, ,  to  B.  should  increase  jointly  with  the  C  capability,  B„  , 

U  d  C 

and  the  number  available  for  change,  Bu-  On  the  other  hand,  (2.1, a)  does  not 

3 

reflect  processing  constraints:  if  each  C  -equivalent  Bc  -  unit  can  service  one 
Bu  -  unit  at  a  time,  then  the  appropriate  function  should  resemble 

cua(Bc(t),  ^(0)  =  euamin^l<Bc^t^  ’  Bu(t^;  (2.1,b) 

here  K  represents  the  conversion  factor  that  allows  combat  units  (e.g.  tanks) 
to  be  interchanged  for  sensor-communication  units,  or  "channels"  for  short.  Other 
forms  for  the  conversion  rate  Cua  may  be  derived,  possibly  by  modeling  this  oper¬ 
ational  segment  in  the  light  of  empirical  study  of  any  data  that  happens  to  be 
available. 

Together,  the  two  equations  (2.1)  and  (2.2)  simply  state  that  a  decrease  in 
B  (t)  translates  into  an  increase  in  B.(t)  during  the  time  period  (t,  t+dt). 

u  d 

For  simplicity,  there  is  no  attrition  of  B  by  R  ,  except  for  that  allowed  to 

3 

deplete  B's  C  capability,  Bc(t). 

The  equation  (2.3)  states  that  the  rate  of  decrease  of  the  Blue  C3  force, 

3 

B  ,  depends  upon  the  magnitude  of  Red  Counter-C  activity,  Rrr(t),  as  well  as 
t  c  c 
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f 
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Blue's  C  activity,  Bc( t ) .  Tentatively  illustrate  by  the  simple  relationship 

Ccc(Rcc(t))  ■  -  ♦•Rcc(t),  (2.3,,) 

<p  being  a  constant.  This  Is  a  conventional  Lanchester  aimed-fire  model.  Another 
similar  model  might  better  reflect  R's  lack  of  knowledge  of  Bc  location: 

Ccc("cc<t»  '  -  ♦  •feE)  '  Rcc(t)  (2.3,b) 

where  BQ  is  proportional  to  the  area  in  which  Bc  is  concentrated,  and  on  which 
Rcc  concentrates  activity.  This  is  essentially  an  unaimed  fire  model  of  classical 
Lanchester  vintage.  There  may  be  reason  to  transition  from  (2.3,b)  to  (2. 3, a) 

3 

during  the  course  of  the  engagement,  as  information  about  Blue  C  increases. 

•  Blue  and  Red  Physical  Attrition. 

TT*  '  “cc(Rcc<'>’  "„<*>•  Ba(t»  <2'4> 

<K» 

"dT  =  ”  °A^A^*  Bu(t),  Ba( t) )  (2.5) 

The  rate  functions  “cc(*)  and  a^(-)  represent  the  physical  attrition  exacted  by 

3 

R  on  the  two  identified  B-force  components:  the  C-C  force  Rcc  ,  and  the  active 
attacking  force  R^  .  It  is  the  magnitude  of  R^(tg),  i.e.  the  attack  force 
survivorship  at  the  time  tg  when  the  bastion  is  reached,  that  determines  whether 
the  bastion  is  actually  taken.  However,  the  size  effectiveness  and  vulnerability 
of  Rcc  can  indirectly  but  decisively  influence  the  latter  variable  and  hence 
the  conflict  outcome.  The  interplay  of  the  variables  described  by  equations 
(2.1)  -  (2.5)  quantifies  the  qualitative  system  behaviour. 
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Here  are  some  specific  expressions  for  physical  attrition 


»cc(Rcc<t>’  Ba(t»  ■ 


/*cc<t)V 

‘“•“YW 


B  (t)  +  D 
u  '  a,cc 


kRcc(t) 


kRcc(t)  +  RA(t) 


Ba(t) 


(2. 4, a) 


The  first  right-hand-side  (rhs)  term  represents  the  rate  of  attrition  of  Rcc  by 
unaimed  B  fire;  the  parameter  p  is  composed  of  both  B's  fire  rate  and 

U  jLL  U 

single-shot  kill  probability.  The  second  rhs  term  represents  the  rate  of 
attrition  of  R__  by  aimed  B  fire;  the  parameter  p  summarizes  the  joint 
effect  of  Bfl  firing  rate  and  kill  probability.  The  ratio  k  Rcc/kRcc  +  RA) 
expresses  the  fraction  of  the  aimed  fire  that  is  directed  at  the  Rcc  force, 
where  k  is  a  parameter  that  may  be  adjusted  to  account  for  various  combat- 
related  effects,  for  example:  the  relative  hardness  or  invulnerability  of  the  Rcc 
force  segment  as  compared  to  the  RA  force  segment.  It  also  accounts  for  the 
relative  exposures  of  the  two  forces  (Rcc  and  RA)  to  aimed  fire  by  Ba  .  Note 
that  if  k  ■  1  the  probability  that  a  unit  of  aimed  fire  will  be  directed  at 
the  Rcc  element  is  simply  RCC/(RCC  +  RA)  —  the  fraction  of  the  existing  P 
force  at  time  t  that  is  devoted  to  C-C  .  If  k  =  0,  no  aimed  shots  are 

3 

directed  at  the  C-C  element,  instead  being  concentrated  on  the  attackers,  RA  , 

3 

while  if  k  -*  °°  B  aimed  fire  is  concentrated  on  Rcc  ,  the  C-C  element.  In 
short,  the  simple  parameter  k  expresses  the  capabilities,  vulnerabilities,  and 
priorities  of  both  B  and  R  .  It,  or  a  more  elaborately  developed  counterpart, 
represents  combat  decision  choice,  and  would  appear  to  have  important  influence 
on  the  progress  of  the  combat. 

The  rate  of  decrease  of  the  attack  force  may  be  specifically  expressed  as 
follows: 
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aA(RA(t).  Bu(t),  Bg(t)  = 


P 


uA* 


RA(t) 

Bu{t)  +paA-  k  ftcc(t)  7  RA(t)_ 


Ba(t) 


(2. 5, a) 


Again  the  first  rhs  term  represents  the  attrition  rate  component  resulting  from 
the  relatively  ineffective  (unaimed)  B  fire.  The  second  rhs  component  repre¬ 
sents  the  attrition  rate  from  aimed  fire  by  B.  The  fraction  RA/(k  Rcc  +  R  A) 
is  simply  the  complement  of  that  appearing  in  (2. 4, a);  it  expresses  in  the  simplest 
way  a  composition  of  the  relative  vulnerabilities  of  R  components  and  the 
priorities  of  B. 

There  appears  to  be  little  hope  of  obtaining  insights  directly  from  the 
differential  equations  (2.1)  -  (2.5).  Of  course  such  equations  can  be  solved 
numerically,  as  was  done  for  the  GT  model  by  Tongue  (1979).  But  another,  some- 

J 

what  simpler,  alternative  is  to  express  the  functions  directly  in  discrete  time, 
taking  the  time  steps  to  be  of  unit  size;  perhaps  0.25  hr.  might  be  appropriate 
for  a  start. 


Lanchester-Style  Equations  in  Discrete  Time 

•  Blue's  and  Physical  Interactions;  Red's  C-C^ 

Bu(t  +  1)  =  Bu(t)  -  Cua(Bc(t),  Bu(t))  (2.6) 

Ba(t  +  1)  =  Ba(t)  +  Cua(Bc(t),  Bu(t))  (2.7) 

Bc(t  +  1)  =  Bc(t)  +  Ccc(Bc(t),  Rcc(t) )  (2.8) 

These  are  easily  seen  to  be  the  counterparts  of  (2.1)  -  (2.3).  They  may  be 

solved  recursively,  starting  with  the  initial  conditions  specified  by 
B, ,  ( 0 ) ,  B,(0),  and  B^(0).  Choice  of  the  initial  condition  by  B  constitutes  a 
decision,  for  if  initial  force  size  is  B  =  100  ,  then  choosing  Bc(0)  =  20 

3 

leaving  B,,(0)  =  80  and  B,(0)  =  0  implies  more  faith  by  B  in  his  C  capabilities 

U  a 

and  invulnerability  than  does  the  choice  B  (0)  =  50  leaving  B,,(0)  =  50  with  B _ ( 0 )  =  0. 

C  U  a 

The  same  considerations  hold  true  for  R  also. 
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•  Red's  Physical  Attrition 

"celt  +  »  -  *m<‘>  -  “ccOccl1)-  Bu(t)-  Ba(t))  (2'9) 

RA(t  t  ))  ■  RA(t)  -  aA(RA(t>,  Bu(t),  Ba(t)  (2. >0) 

Here  again  one  starts  with  initial  conditions  Rcc(0)  and  ( 0)  and  solves 

recursively  to  find  the  Red  attacker  force  size  at  tQ,  when  the  defended  position 

3 

is  reached.  If  R  has  considerable  faith  in  the  unit  effectiveness  of  its  C-C  , 
then  presumably  R  (0)  is  chosen  to  be  relatively  small,  permitting  the  majority 
of  its  resources  to  oe  allocated  to  the  attack  force  R^(0). 
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3.  NUMERICAL  ILLUSTRATIONS  OF  MODEL  I  PERFORMANCE 


The  present  version  of  Model  I  is  as  simple  as  seems  consistent  with  our 
attempts  to  blend  elements  of  Conwand  and  Control  with  combat  interactions 
between  forces.  Even  so  five  state  variables  are  needed  to  describe  system  be¬ 
havior,  and  their  inter-related  evolution  in  time  is  sufficiently  complex  to 
make  a  direct  mathematical  discussion  appear  unprofitable.  As  an  alternative 
we  have  elected  to  create  a  computer  program  that  produces  the  numerical  se¬ 
quences  of  values  assumed  by  the  various  forces  as  combat  progresses. 

The  computer  program  is  written  in  FORTRAN;  a  listing  appears  in  an  Appendix. 

If  combat  is  terminated  at  a  particular  time  point,  tg,  the  relation¬ 
ship  of  R^(t)  to  B^(t)  will  be  assumed  to  determine  the  outcome.  The 
latter  relationship  is  itself  influenced  by  decision  variables  on  each  side. 
Here  are  some  options  and  constraints  for  the  combatants. 

•  Initial  Red  attacker  force  size,  R(0)  ,  and  its  division 
into  attackers,  R.(0);  and  Counter-C^  (C-C^),  R  (0);  then 

A  C  C 

R(0)  =  RA(0)  +  Rcc(0) 

•  Initial  Blue  defender  force  size,  B(Q),  and  its  disposition 
into  Uneffecti ves/Undesignated,  B  (0);  Effectives,  B  (0); 

U  a 

and  C3  Forces,  B  { 0) ; 

B ( 0 )  =  11,(0)  +  B  (0)  +  8  (0) 

U  d  C 

Often,  B  (0)  :  0  will  be  reasonable  as  an  initial  condition;  this  would 
a 

represent  surprise  by  the  Red  force. 

3 

Note  that  Blue  has  a  technological  constraint  that  limits  its  C 

capacity.  Two  parameters  actually  play  this  role:  (1)  0,  .  representing 

the  rate  per  unit  time  of  acquisition  and  transfer,  and  thus  of  converting 

B  's  into  B  's,  and  (ii)  K,  the  capacity  factor  according  to  which  force 
u  a 

3 

units  ( B  ( t )  =  tank  equivalents,  say)  are  made  equivalent  to  C  units 
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(sensor-communication  combinations  or  "channels"  for  short);  see  (2.1,b)  for 
the  rate  expression  actually  used.  Initially  imagine  Blue  to  be  merely  en¬ 
dowed  with  these  parameters,  and  allow  them  to  remain  fixed.  It  may  be  rea¬ 
sonable  for  them  to  change  as  combat  progresses  and  sensors  are  disabled, 
etc.  The  current  computer  program  can  be  straightforwardly  altered  to  reflect 
such  combat  related  damage  or  degradation. 

Blue  is  assumed  to  be  in  possession  of  the  decision  parameter,  k, 
see  (2. 4, a)  and  (2. 5, a):  increases  in  k  directs  a  greater  proportion  of 
B  's  fire  at  R  forces,  while  reduction  of  k  concentrates  B,  fire  at 
the  Red  attacking  force,  R^;  k  =  0  means  exclusive  concentration  on  RA  . 
Although  the  present  model  explicitly  makes  k  a  constant  throughout  the 
combat  period,  no  such  restriction  need  be  at  all  permanent:  it  may  perhaps 
be  best  for  B  to  switch  from  k  =  °°  at  early  stages  of  combat  --  thus 
maximally  reducing  Rcc's  interference  with  initial  rapid  buildup  of  the 
Ba‘s  force  —  finally  switching  to  k  =  0  later  on  so  as  to  concentrate  on 
decimating  the  R^  force  before  the  defended  position  is  reached.  If  for 
example,  Bc(t)  is  (i)  hardly  reduced  at  all  after  a  few  periods,  or 
(ii)  is  almost  wiped  out,  then  there  would  seem  to  be  little  reason  for  Ba 

a 

tc  target  Rcc  units  any  longer.  Of  course  if  Bu  forces  are  almost 
entirely  converted  to  B  's  there  would  again  seem  to  be  little  reason  to 

a 

target  Rcc>  for  Bc's  function  has  been  accomplished  --  at  least  so  far  as 
the  current  engagement  is  concerned. 

Switching  betwen  two  extreme  values  of  k  is  reminiscent  of  the  "bang- 
bang"  policies  of  control  "optimal,"  but  there  is  no  measure  of  effectiveness 
or  figure  of  merit  yet  specified  for  B. 
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Numerical  Cases 


Here  are  the  parameter  values  selected  for  initial  exercise  of  the  model 
Ra(0)  =  120  >  Rcc(0)  =  30 

<t>,  (attrition  rate  of  Bc  by  R  )  ,  =  1  . 

Ba(0)  =  0  ,  Bu(0)  =  80  ,  Bc(0)  =  20 

puA  ”  pu,cc  ”  ® 

For  simplicity  we  are  assuming  that  the  Unaffective  Blues  have  no  combat  per¬ 
formance  capability,  but  are  merely  a  pool  of  assets  from  which  Effectives  are 

3 

created  by  the  information  from  the  central  C  facility.  This  assumption 
can  easily  be  modified  if  desired. 

Note  carefully  that  in  all  that  follows  the  numerical  values  for  various 
parameters  have  been  chosen  for  illustrative  purposes  only,  and  need  bear  no 
close  resemblance  to  any  actual  values,  which  are  in  fact  unknown  to  us.  The 
purposes  of  the  cases  discussed  is  entirely  exploratory.  But  of  course  we  hope 
that  that  the  suggestions  and  implications  noted  will  promote  interest  in 
further  work,  leading  to  model  refinements  and  the  use  of  more  nearly  correct 
parameter  values,  if  such  can  be  agreed  upon. 
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C*se-J :  Pa,cc  =  1-° 
$  ■  1.0 


paA  *  1*0  (Blue  Attrition  rates) 

3 

(Red  Attrition  rate  on  Blue  C  ) 


Medium  Speed  Blue  C3 

Target:  Attack  C-C3 

Rapid 

Blue  C3 

Attack 

C-C3 

m 

11 

CM 
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II 

k=0.2 

k=5 

eua~^  9ua'"^ 

0  =2 
ua 

eua"‘ 

It 
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0.00 


3.0 

1.50 

0.60 

0.00 


6  0.00  .06 

7 

0.00 

8  ' 

*  4^  > 

Red/Blue  Active  Force  Ratio  as  Combat  Progresses 
(Red  C-C3  Attrition  Rate  <p  =  1.0) 

Fig.  3.1 


The  numbers  in  the  above  figure  suggest  that 

3 

t  Increasing  the  speed  e  ,  (reducing  the  response  time)  of  Blue  C 
activity  has  a  profound  effect  upon  the  Red/Blue  Force  Ratio: 
e.g.  at  t  -  4  when  k  *  0  (Blues  concentrated  on  Red  3 
Attackers)  Reds  outnumber  Blues  by  1.6  to  1  if  Blue’s  CJ 
speed  is  e  *  1,  while  doubling  that  rate  to  eya  *  2  under 

otherwise  the  same  circumstances  cuts  the  ratio  to  a  Blue- 
favorable  0.21  to  1,  l.e.  by  a  factor  of  about  eight. 

•  For  the  parameters  considered,  concentration  by  Blue  on  Red 
Attack  forces  pays  off  more  than  does  concentration  on  the 
Red  C-C3:  e.g.  at  t  =  4  with  Rapid  Blue  C3  the  advantage  to 
Red  changes  from  a  force  ratio  of  1.6  to  2.1  with  an  increased 
concentration  of  fire  by  Blue  on  Red's  C-C3. 
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The  comparative  advantage  of  Red  to  Blue  (RA(t)/Ba(t) ,  yielded 

by  Blue  emphasis  on  Red  Attackers  (k  =  0.2)  instead  of  Blue 

emphasis  on  Red  C-C  (k  =  5)  diminishes  as  Blue  C  speed  e  , 

decreases.  The  reason  may  be  that  when  the  crucial 

Blue  function  is  relatively  weak  or  slow,  it  is  profitable 
to  spend  more  effort  n  its  defense,  at  least  until  it  has 
served  its  assigned  purpose  of  converting  B's  to  B ,'s  . 

Case  2:  Pa>cc(t)  =  0.25  +  0.075t  =  pflA(t) 

=  1.0 

In  this  case  a  space  (and  hence  time)  varying  representation  of  the  attrition 
rate  is  introduced:  at  t  =  0  when  the  Red  force  is  far  away,  attrition 
rates  on  it  are  taken  to  be  low,  but  they  increase  steadily  with  decreasing 
range  (increasing  t). 


Medium-Speed  Blue  Cv 
Attack  C-C^ 

k=0. 2  k=5 


Rapid  Blue  C" 


Attack 
k =0.2 


It 

eua=1 

w1-0 

6  =2.0 
ua 

e  =2.i 
ua 

1 

03 

00 

CO 

oo 

2 

6.0 

6.00 

3.00 

3.00 

3 

3.35 

3.44 

1.58 

1.60 

4 

2.25 

2.50 

0.88 

1.12 

5 

1.37 

1.82 

0.23 

0.60 

6 

0.40 

0.89 

0.00 

0.00 

7 

0.00 

0.00 

si 

si 

Red/Blue  Active  Force  Ratio  as  Combat  Progresses 
(Red  C-C^  Attrition  Rate  <j>  =  1.0) 

Fig.  3.2 
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There  is  little  qualitative  change  in  the  numbers  obtained,  as  compared  to 
those  of  Fig.  3.1.  The  change  in  the  attrition  rate  function,  increasing  from 
a  small  to  a  larger  value  as  time  goes  on,  allows  a  larger  number  of  Reds  to 
survive  longer  than  was  true  for  Case  1.  Again  it  seems  profitable  for  Blue 
to  expend  most  of  its  energy  on  Red  Attackers. 

"a,cc  =  °-25’  -aA  =  °-25 

<P  =  1.0  (Red  Attrition  Rate  vs_  Blue  C2) 

(■•  =  O.'i  (Relatively  Slow  Blue  C^) 

j  Jd 

Tua  =  (Relatively  Rapid  Blue  C2) 

In  order  to  describe  the  results  the  entire  computer  printout  is  now  presented. 
Parameter  values  are  shown  across  the  oage  top;  note  that  KC  (computer 
printout)  =  k  (text);  the  larger  (smaller)  this  parameter  becomes  the  greater 
(less)  is  concentration  of  Elue  fire  on  Red  C-C  . 

•  Compare  cases  for  which  r>  =  0.5  to  0113  -  1.0:  the  force 

U  a  U  a 

nation  R./B,  for  the  slower  system  (0iia  =  0.5)  is  about  twice 

Ma  lid 

that  for  the  faster  (t>  =  1.0)  system,  no  matter  which 

ua 

faring  strategy  is  adopted  (k  =  0.2.  5,  or  10). 

•  The  time  period  during  which  R^/ it  above  unity  (Reds 

have  advantage),  or  there  is  near  parity,  is  about  15  to  16 

3 

for  6ua  =  0.5  iRelaMvely  Slow  Blue  C  );  the  same  period  is 
about  8  to  9  for  =  1.00.  After  that  period  the  R. 

Ud  H 

force  is  quickly  wiped  out,  but  there  are  wide  differences 
in  the  readiness  of  tne  8,'s  :  for  k  =  KC  =  5,  10  the  B  's 

a  d 

soon  reach  their  maximum  of  80,  (t  ~  12),  while  at 
k  =  KC  =  0.2  the  maximum  is  not  reached  until  t  ~  42. 


^  •*.  04  U4040404CUC40404  04  MM  M  !SJK>N>rO  Mf4M  *■*  ►**-*►*  *+  *-►***►-,  t— 
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Red/Blue  Active  Force  Ratio  as  Combat  Progresses 

Fig.  3.4 

The  most  striking  effect  visible  in  the  table  is  that  Blue  improves  his 
performance  relative  to  Red  by  decreasing  the  Effectives  and  increasing 
the  allocation  of  his  forces  to  C  .  A  nearly  3-to-l  improvement  (for  B)  of  the 
force  ratio  at  t  =  5  is  apparent.  The  effect  of  the  parameter  k,  which 
dictates  the  fraction  of  B  energy  expended  to  deplete  Red  C-C  ,  is  very  small. 
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4.  A  STOCHASTIC  VERSION  OF  MODEL  I. 

The  evolutionary  equations  (2.1)  -  (2.5),  or  equivalently  (2.6)  -  (2.10) 
are  entirely  deterministic,  a  feature  that  seems  unrealistic  since  in  reality 
uncertainty  and  random  variability  abound.  There  are  several  ways  in  which 
uncertainty  may  be  allowed  to  intrude  into  the  formulations;  for  instance 

(i)  Through  the  necessity  of  estimating  parameters 

(plt  rr,  p,  k,  etc.)  from  data,  assuming  the  model  specification 

UfCC  a  ) C 

is  "correct",  or  at  least  adequate. 

(ii)  Through  the  necessity  of  using  simplistic  models  such  as 
(2.1)  -  (2.5)  to  represent  a  more  complex  reality. 

(iii)  By  explicitly  permitting  randomness  to  enter  the  dynamic  equations 
as  an  additional  driving  force. 

(iv)  Other  possibilities;  combinations  of  the  above,  for  example. 

Begin  by  adding  a  random  perturbation  term  to  the  discrete  time  equations 
(2.6)  -  (2.10)  as  in  (iii)  above: 

vt+,>  ■  v11  -  MV11,  v*>)  ♦  °uK*i(u>) 

4a<»+'>  *  Vl>  *  Cua(Vt>'  V*>)  *  0aK*1(a>) 

*c<t+l>  *«c(t>  -  C^U),  Rcc(t)  ♦  «c  4Wtt,(c>)  (4.1) 

4c<t+,>  '  &cc(t)  -  “cc(«cc(t>-  V1)-  &,<«>)  +  0cc(a*Wcc)) 

JBfttt+l )  -  %(t)  -  aA(Rfl(t),  ^(t),  6a(t))  +  ofl(iKttl (A)) 

The  notation  j(  denotes  a  random  function.  The  vector  of  random  components 
j^(u),  J^(a),  ^(c),  ^(cc),  ^(A)  will  be  taken  to  be  one  of  not  necessarily 

independent  Wiener  processes,  sampled  at  time  points  t  =  1,  2,  3,  -  Thus 

the  equations  (4.1  )  turn  out  to  be  analogous  to  Ito-type  stochastic  differential 
equations;  such  equations  have  been  used  by  Lehoczky  and  Perla  (  1978)  to  describe 
combat  situations. 
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J 


¥ 


4.1.  Explicit  Representation  of  Random  Terms. 


By  arguments  analogous  to  those  of  Lehoczky  and  Pert  a  (  1978  ),  or  of  Gaver 
and  Lehoczky  (  1977)  in  a  different  context,  we  write  down  expressions  for  the 
scales  o ( • )  of  the  random  components  written  as  o.  ^AWt+i  ( •  )j  (4.1). 

We  argue  heuristically  that  if  the  stochastic  processes  are  nearly  Markovian  and 
further  are  superpositions  of  many  point  (e.g.  birth-death)  processes  describing 
the  changes  of  individuals  states  in  a  relatively  short  time  period,  then  it  is 
reasonable  that  o ( • )  be  equal  to  the  square-root  of  the  individual  drift 
(deterministic)  terms,  the  latter  being  given  by  the  expression 

Cua(V^’  V^)’  CccM>*  V^)’  etc-  aPPearin9  on  the  rhs  of  (4.1). 

The  latter  is  suggested  by  the  approximately  Poisson  nature  of  the  changes 

in  state  of  the  system  over  relatively  short  time  intervals.  Recall  that  the 
standard  deviation  of  a  Poisson  random  variable  (state  change,  here)  is  equal 

to  the  square-root  of  the  mean  (state  change),  or  drift  in  diffusion  theory  jargon. 
Furthermore,  the  random  elements  AWt+^(n),  AWt+^(a),  etc.,  are  realizations  of 

independent  unit  normal  or  Gaussian  random  variables.  Some  readers  will  recognize 
tne  resulting  system  to  be  a  discrete-time  version  of  Ito-type  stochastic 
differential  equations;  see  Arnold  (  )  for  basic  information. 

Mere  are  some  explicit  examples,  following  (2.1,b),  (2.3),  (2. 4, a),  (2. 5, a). 


wnere 

AWt+-|(a)  ~N(G,1),  i.e.  has  the  Normal  distribution  with  mean  0, 

and  standard  deviation  1, 


ana 

AW^-j(a)  -  -  AW^-jiu), 
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the  latter  because  a  fluctuation  away  from  the  mean  change  in  one  di recti on--down, 
say--for  B  (t)  is  exactly  matched  by  one  in  the  up  direction  for  B  (t)  during  time 

U  a 

interval  t  to  t+1 .  Next,  using  (2.3,b), 


v(Awt+i(c)) = b0  )Rcc^t)  (Awt+i(c0 


(4.3) 


where  AWt+^(c)  ~N(0,1),  and  AWt+^(c)  is  independent  of  AWt+^(a)  as  well  as 

of  past  values  aWt(c),  t  =  0,1,2 . t. 

Next , 


u,cc 


LR0  3 


Bu(t) 


(AWt+1(u,Cc)) 


(4.4) 


/ ' 

■  KRcc(t)  1 

D  /  +  \ 

a ,cc 

kRcc(t)  7  RAlt)_ 

Ba{t) 

(AWt+l (a*cc)) 


where  AWt+-j(cc,u)  and  AWt+^(cc,a)  are  independent  random  elements,  both~N(0,l), 
that  represent  respectively  the  fluctuation  away  from  mean  attrition  on  the  Red  C-C 
facility  caused  by  unaimed  (ineffective)  fire  and  by  aimed  (effective)  fire. 
Finally,  use  of  (2.5a)  provides 


vK+i<a>) 


|W| 

B  (t) 
uv  ' 

duA 

R0 

(a*WU*A)) 


(4.5) 


1 

— 1 

VPaA‘ 

nrxtrn^tT 

Ba(t) 

- 

^AJd^-]  (a,A)j  . 


Again,  aW^.+^(u,A)  and  aW^-jU.A)  are  independent  and  normally  distributed,  rep¬ 
resenting  fluctuations  around  mean  attrition  on  the  -  component.  Under  the 
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present  simple  circumstances  in  which  approximate  independence  may  be  justified, 
we  can  actually  combine  the  two  "noise"  terms  in  (3.4)  and  (3.5)  to  obtain  the 
simpler  forms 


~cc’  'J^'t+l^cc' 


u,cc 


‘Rcc^n 


0  J 


B  (t)  +  p 
it  '  a,cc 


kR  (t) 
ccv  ' 


kRcc(t)  +  RA(t). 


Ba(t) 


(-i<it+i (cc) 


(4. 4, a) 


zk  ^+1 


(A) 


rr 

j 

.  / 

\w] 

V'uA 

[R0  J 

Bu^  +  paA 


RA(t) 


kR"cc(t)  + 


Ba(t) 


+  1 


(A) 


(4.5, b) 


s  will  reduce  simulation  difficulties  by  making  it  unnecessary  to  sample  two 
normal  random  variables  for  each  of  the  last  equations.  The  introduction  of 
correlation  is  easy-  if  justified. 

It  might  be  pointed  out  that  considerable  freedom  exists  in  the  choice  of 
tne  distribution  of  all  of  the  noise  terms  aW^.+-.  ( • ) :  they  need  not  be  normal, 
•.o**  need  tney  be  independent,  nor,  in  fact,  need  they  be  independent  of  the  cor¬ 
responding  state  value.  Of  course  a  random  fluctuation  that  sends,  say,  JR,^ (17) 

*■■  ■  a . '/  regat'.e  must  in  truth  merely  wipe  out  the  A-force,  i.e.  reduce  it  to 
Intuition  indicates  that  a  formal  (and  forbidden)  passage  below  zero  for 
;;  1 t-I),  and  that  the  latter  is  unlikely  to  be  extensive  by  virtue  of  the  small 
rO‘se  variation  near  £ ,  =  0  for  the  particular  form  of  (4.5,b).  However,  choice 
of  a  more  gaudily  variable  noise  increment,  e.g.  with  AW. +, ( • )  now  chosen  to  be 
•w.  g  tailed,  perhaps  in  a  Cauchy-like  manner,  will  likely  require  more  extensive 
fixing  at  the  boundary.  This  is  no  reason  to  avoid  such  models,  for  there  is 
•  ctnir.c  *' o ’ y  aoout  tre  normal  (or  Wiener-like)  variation  save  for  its  appearance 
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as  the  noise  when  a  basic  Markov  structure  is  assumed  to  underlie  the  present 
models.  Actually,  over-variation  (from  the  Gaussian/normal)  and  serial 
dependence  may  well  usefully  represent  mixtures  of  normals  (or  other)  dis¬ 
tributions  resulting  from  factors  such  as  terrain,  visibility,  and  many  other 
features  which  combine  to  generate  departures  from  the  systematic  determin¬ 
istic  models  analogous  to  the  classical  Lanchester  forms  (2.1)  -  (2.5). 
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5.  MODEL  II:  TWO  FORCES  IN  COMBAT  WITH  MUTUAL  ATTRITION 


The  model  described  earlier  dealt  with  a  special  situation  in  which  one  force 

(Red)  was  the  attacker  of  another  in  a  defended  position  (Blue).  Lack  of  symmetry 

was  evident.  The  model  of  this  section  represents  an  equally  stylized  but  now 

symmetric  situation  in  which  both  forces  are  capable  of  causing  attrition  on  each 

other.  Once  again,  guidance  is  furnished  by  explicitly  represented,  and  vulnerable, 

3 

C  agencies.  Our  models  allow  for  different  targetting  strategies. 

Notation  is  as  follows: 

R  (t)  =  the  number  of  acti ve  Reds  at  t,  meaning  the  number  of  R's 
actually  targetted  on,  and  firing  at,  Blue  units. 

Bit)  =  the  corresponding  Blue  force  size. 

a 

R. (t)  =  the  number  of  ineffective  (or  inacti ve)  R's  at  t,  meaning 
1  the  number  of  Rrs  currently  untargetted  and  awaiting  new 

assignment. 

B. (t)  =  the  corresponding  Blue  force  size. 

3 

R  (t)  =  the  size  of  the  Red  C  agency  at  t,  i.e.  the  force  responsible 
c  for  switching  R.  -  units  to  Ra  -  units. 

Bc(t)  =  the  corresponding  Blue  force  size. 

The  quantities  R  (t)  and  6c(t)  are  to  be  viewed  as  the  command  and  control 
authorities  responsible  for  the  individual  direct  combat  elements  on  their  respective 
sides.  The  state  variables  R,(t),  R.(t),  B  (t),  B.(t)  may  be  thought  of  in 

did! 

units  of  individual  tanks  or  ships,  or  as  aggregations  such  as  battalions,  companies, 
naval  task  groups,  or  whatever  is  appropriate  for  the  particular  situation  under 
consideration.  The  intention  of  the  present  model  is  to  simply  express  changes 
in  the  respective  force  sizes  in  terms  of  rate  processes--both  information  transfer 
rates  and  physical  attrition  rates--and  in  terms  of  initial  allocations  of  resources. 

It  is  clear  that  when  Red  and  Blue  forces  come  into  contact  a  variety  of  pos¬ 
sible  behaviors  may  occur.  The  models  developed  here  are  intended  to  represent  the 


24 


consequences  of  a  few  simplified  versions  of  the  true  complexity  possible.  In 
particular,  they  permit  the  study  of  different  target  category  priority  schemes. 

Red(Blue)  Force  State  Equations 

The  general  form  of  the  state  change  equations  now  follows.  They  are 
given  only  for  Red,  but  the  Blue  equations  are  symmetric. 

Ra(t+1 )  =  R*(t)  -  CBRA(R*(t),  R*(t) ,  R*(t) ,  B*(t) ) 

^  ^  ^  -  - —  ■  S 

Physical  Attrition,  R.  (5.1) 

a 

+  DRA(R*(t),  R*(t) ,  R*(t),  B*(t),  B*(t),  B*(t) ) 

Effect  of  Information  State  Change,  Ra 

R1 (t+1 )  =  R*(t)  -  CBRI(R*(t),  R*(t) ,  R*(t) ,  B*(t)) 

Physical  Attrition,  R^  (5.2) 

+  DRI (R*( t) ,  R*(t) ,  R*( t) ,  B*(t),  B*(t) ,  B* ( t ) ) 

v - - - - , - ' 

Effect  of  Information  State  Change,  R^ 

Rc(t+1)  =  R*(t)  -  CBRC(R*(t),  R* ( t ) ,  R*(t),  B*(t) ) 

Red  C  Capacity  Attrition  and  Suppression  (5.3) 

+  DRC(R*(t),  R*(t) ,  R*(t),  B*(t) ) 

_ _ _ y  ■  ■  - - ' 

3 

Red  C  Capacity  Restoration  and  Recovery 

★  ★  ★ 

\  •  Additionally,  to  compute  Ra,  R. ,  R.  from  R  ,  R . ,  R  ,  utilize  the  appro- 

1  a  I  C  a  1  C 

priate  physical  constraints: 
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R*(t+1)  =  min[Ra(0)  +  R.(0) 

R*(t+1 )  =  min^Ra(0)  +  R^O) 
R*(t+1)  =  mi n |R  (0),  max{R 


max{R,(t+l ) , 

a 


0} 


,  max{R..  (t+1 ) ,  0}J 

(t+1),  0}] 


•  Initial  conditions,  i.e.  values  of  R  (0),  R,.(0),  R  (0),  are  required  to 

a  1  C 

start  the  process,  after  which  the  various  force  sizes  are  calculated 
recursively  in  time. 

•  Abbreviate:  CBRA(arguments  as  in  (5.1))  s  CBRA(t) 

[^(arguments  as  in  (5.1))  =  0RA(t), 


etc. 

Blue  State  Equations  resemble  (5.1),  (5.2),  (5.3)  with  changes  in  rate 
definition  and  notation 

Red  Blue 

CBRA  CRBA^Ba^’ 

DRA  dba 

CBRI  CRBI 

DRI  dbi 

CBRC  CRBC 

dbrc  drbc 


There  now  follow  some  specific  expression  for  the  above  rates,  intended  to 
illustrate  possible  effort  allocation  strategies. 

1 


Hr 

■f 
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A.  Greedy  Allocation  (Blue  vs.  Red) 

3 

Priority  Order:  Red  Actives,  Red  Inactives,  Red  C 

a)  CBRA^  =  PA  * 

B ( R_ ( t ) )  =  Number  of  Blue  Actives  vs.  Red  Actives 

C 

=  min^B*(t ) ,  Ra ( t )J  (5.4) 

•  Note  that  all  B  's  fire  one-on-one  onlR  's  if 

-  a  a 

B,  force  <  R,  force;  otherwise  one-on-one  until  R,  targets 

a  a  a 

are  Insufficient,  leaving  BJt)  -  RJt)  to  be  used  against 

a  a 

R^'s,  which  are  next  on  the  priority  list. 

b)  CBRI(t)  =  PI 

★ 

B ( Rj ( t ) )  =  Number  of  Blue  Actives  vs^.  Red  Inactives 

=  minjmax{Ba(t)  -  B(Ra(t)),  0},  (t (5.5) 

C)  CBRC(t)  =  pc  B(R*(t) ) 

B(R*(t))  =  Number  of  Blue  Actives  vs.  Red  C3 

=  min[max{Ba(t)  -  B(R*(t))  -  B(R*(t)),  0),  R*(t)]  (5.6) 

•  Change  of  priority  is  easily  accomplished:  if  priority 
sequence  is 

x  first,  y  second,  z  third 
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cBRx(t)  =  px  B ( R* ( t ) ) ,  B ( R* ( t ) )  =  m1n|B*(t) ,  R*(t) 


CBRy( t )  =  py  B ( Ry ( t ) ) ,  B(R*(t) )  =  min  max{B*(t)  -  B(R*(t)),  0}, 

4.  n 


R*  ( c ) 
y 


5.7) 


CBRz(t)  =  pz  B(Rz(t))’  B ( R* ( t ) )  =  min 


R*(t) 


max)Ba(t)  -  B(Rx(t))  -  B(Ry(t)),  0}, 


•  Same,  with  obvious  changes,  for  Red  on  Blue 

•  Attrition  rates  of  Blue,  by  Red,  are  denoted  respectively  by  6,,  6,,  6.. 

•  Attrition  rates  can  be  made  time  (space)  dependent. 

•  Note  that  the  attrition  law  is  assumed  for  illustration  only  to  be 

in  accordance  with  a  Lanchesterian  linear  law.  Other  types  of  laws 
can  clearly  replace  this  one. 


Proportional  Allocation  (Blue  vy.  Red) 

kA  <U) 


a) 


CBRA(t)  =  PA 


kA  Ra(t)  +  kj  R1  (t)  +  kc  Rc(t) 


B>) 


(5.8) 


•  k^,  kj.  k^  are  control  "constants",  representing  target  priorities. 
They  can  be  changed  or  combat  progresses,  if  desired.  If  k^  |  °° 
(practically  k^/k j  ~  k^/k^  ~  10  should  do)  there  is  heavy  emphasis 
on  R  1  s 

d 

•  is  attrition  rate  (firing  rate  times  single-shot  kill  probability). 


b) 


CBRI (t )  =  PI 


R*(t) 


kA  Ra(t)  +  kj  R.(t)  +  kc  Rc(t) 


B^tt) 


(5.9) 
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c)  CBRC(t)  *  oc 


kc  Rc(t) 


kA  Ra(t)  +  kI  “i(t)  *  kC  Rc(t) 


B*(t) 


•  Same  for  Red  on  Blue  with  appropriate  notational  changes. 

•  This  is  an  alternative  strategy  that  allocates  some  B  -effort  to 

a 

all  components  of  the  Red  force.  It  samples  rather  than  gulps. 
There  are  many  possible  alternatives. 


Red(Blue)  Information  State  Change  Equations 

Here  are  some  sample  equations  that  are  tentatively  proposed  to  model 

2 

the  effects  of  philosophy  and  capability  of  the  C  component. 

(1)  Tight  Central  Control (Reds) ;  C^  Attrition 

Only  R  's  may  fire  at  Blues.  When  a  Red  firing  engagement 

a 

terminates,  it  returns  to  the  R^  status  (becomes  inactive).  The 
R„  force  acts  to  change  or  return  R. 's  to  R  's.  Thus 

C  la 

(a)  DRA(t)  =  min(CR  KR  R^(t),  CR  R.(t),  R^t)} 

v _ —  .  - _ _ _ — ' 

P 

(R.  ->  R  rate;  C  process) 

1  a 

“jj'RBA^  +  CRBI^  +  CRBC^ 

' - - - - - - - - - - - ' 

(Number  of  Red  vs  Blue  Engagements  Terminated 

in  t^1  Period,  or  Number  of  R  ‘s  Released  to 

a 

Rr) 

(b)  DRI (t )  =  -  DRA(t) 


(5.10) 


(5.11) 


(5.12) 
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(5.13) 


2 

(c)  DR{.(t)  =  0.  (No  Red  C  recovery  from  suppression.) 

•  Note  that  in  this  model  both  physical  attrition  (which  is  permanent) 

2 

and  lack  of  information  (temporary,  unless  C  is  inadequate)  reduce 
Red  effectiveness. 

2 

(2)  Loose  Central  Control;  C  Attrition 

Same  as  (1),  but  retain  only  **»•«:*  term 

(a)  DRA(t)  =  min{CR  KR  R*(t),  CR  R.(t),  R.(t)}  (5.14) 

(b)  Dri ( t )  =  -  DRA(t)  (5.15) 

(c)  DRC(t)  =  0  (5.16) 

•  Note  that  this  model  assumes  that  once  an  changes  to  Ra 
status  it  remains  so  until  physically  removed. 

•  A  model  that  lies  between  (1)  and  (2)  may  be  attractive,  in  which 
case  one  might  multiply  the  second  term  in  (5.11)  by  a  time  (distance) 
dependent  factor  less  than  unity  to  represent  the  fraction  of  the 

2 

terminating  Red  engagements  that  must  have  recourse  to  central  C 
re-direction. 
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6.  NUMERICAL  ILLUSTRATIONS  OF  MODEL  II  PERFORMANCE. 

In  this  section  we  present  some  examples  of  Model  II  behavior.  Once  again 

these  numerical  results  are  obtained  from  a  computer  program  written  in  FORTRAN; 

a  listing  appears  in  an  Appendix.  The  intention  is  to  provide  some  feeling  for 

the  sensitivity  of  the  model  to  changes  in  parameters  and  combat  philosophy  or 

priorities.  Careful  study  of  the  outcomes  will  be  helpful  in  leading  one  to 

comprehend  the  complex  dynamics  of  this  model  and  others  evolving  from  it. 

Discussion.  For  convenience,  consider  the  columns  in  the  Figure  labelled  1  9 

from  the  left;  e.g.  Col.  3  is  headed  Cp=l,(2),  under  pc=0.1,  meaning  that 
2 

Red's  C  rate  is  unity,  the  attrition  rate  of  Blue  on  Red  is  pc=0.1,  and  (1) 
Strong  Central  Control  is  in  use  by  both  antagonists. 

2 

•  Each  combatant's  Active  component  first  grows,  as  C  changes  Inactive 
to  Actives,  and  then  dwindles  by  attrition.  This  is  to  be  expected. 

2 

•  Compute  Col.  2  and  Col.  3:  Strong  Central  Control  tends  to  overload  C 

to  the  detriment  of  Red  Actives.  See  also  Col.  6  and  Col.  7:  same  lesson 

for  Blue  Actives.  A  tradeoff  would  occur  if  S.C.C.  by  Red  was  combined 

with  higher  attrition  rates  for  Red  than  Blue,  Blue  using  L.C.C. 

2 

•  Compare  Col.  2  and  Col.  4:  increasing  (doubling)  C  processing  rate 
under  S.C.C.  provides  a  decided  initial  improvement  for  Red.  but  this 
disappears  later;  the  number  of  Inactives  becomes  small,  so  there  are  no 
resources  from  which  to  draw  to  increase  the  complement  of  Actives. 
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Pc-0.1  Pc=0.5  PC=0.1  PC=0.5 


t 

cR=l  ,(1) 

cR=l , (2) 

crx2,(1 ) 

cr=2,(2) 

cR=l,(l) 

cR=l * ( 2 ) 

cr=2,(1 ) 

cr=2,(2) 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

20 

20 

40 

40 

10 

10 

10 

10 

2 

37 

39 

75 

79 

18 

19 

18 

19 

3 

52 

57 

90 

97 

24 

27 

24 

27 

4 

75 

91 

83 

91 

34 

41 

34 

41 

6 

81 

87 

81 

87 

37 

47 

37 

47 

7 

77 

82 

77 

82 

40 

52 

40 

52 

8 

73 

77 

73 

77 

42 

57 

42 

57 

9 

69 

71 

69 

71 

43 

61 

44 

61 

10 

65 

65 

65 

65 

45 

65 

45 

65 

15 

44 

38 

44 

38 

48 

64 

48 

52 

20 

26 

23 

26 

23 

44 

48 
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36 

25 

16 

13 

16 

13 

35 

38 

21 

26 

Col : 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

Fig.  6.1 

cB  «  1.0,  c  =  0.10 

2 

Priorities:  Actives  -  Inactives  -  C 


(1) :  Strong  Central  Control 

(2) :  Loose  Central  Control 
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Discussion.  This  Fig.  6.2  illustrates  the  effect  of  changes  in  targe*  priority 


under  different  philosophies  of  central  control. 

•  All  else  being  equal,  S.C.C.,  indicated  by  (1),  is  a  handicap  as  illus¬ 
trated  by  comparison  of  Col.  2  and  4  for  Reds,  and  corresponding  Cols.  6 

and  8  for  Blues.  The  penalty  for  Blue  is  higher  possibly  because  of  its 

2 

original  smaller  complement  of  C  assets,  and  hence  its  poorer  traffic 
handling  capability. 

2 

•  Suppose  Blue  reverses  firing  priority,  targetting  C  first,  then  R  , 

3 

finally  R^.  Compare  Cols.  (3)  and  (7)  to  Cols.  (2)  and  (6).  Note  that 

2 

initially  R,  actually  improves  when  C  is  first  priority  because  R, 

a  a 

is  only  lightly  diminished  and  Rc  remains  temporarily  adequate.  It 

2 

appears  that  this  strategy  change  by  B  is  ineffective  because  Red's  C 
facilities  are  ample  enough  to  withstand  the  attack. 

2 

•  Suppose  Red  reverses  firing  strategy,  targetting  Blue  C  first,  then  B  . 

a 

and  finally  R^.  Adopt  loose  control,  (2),  for  illustration.  Note  that 

2 

the  effect  on  Ba  of  Reds  C  -first  strategy  is  initially  small,  but  as 

a 

the  combat  proceeds  B.  force  size  is  considerably  reduced. 

a 

2 

The  last  two  comparisons  suggest  that  primary  attach  on  the  opponent's  C 

force  is  advantageous  when  that  force  is  meager  (or  especially  vulnerable),  while 

2 

such  an  attack  is  actually  counterproductive  in  case  the  C  force  is  adequate. 

Of  course  this  simply  suggests  that  high  priority  is  best  placed  on  attacking 
the  weakest  point  in  the  system.  In  the  present  example  weakness  is  merely  a 

matter  of  numbers  (C^  force  size),  whereas  in  reality  there  must  be  an  assessment  ! 

i 

of  a  potential  target's  capability  or  performance  rate.  Models  of  the  present  typi 
and  their  offspring  and  siblings,  should  be  of  use  for  evaluating  different  pro¬ 
posed  targetting  strategies  or  doctrine. 
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Ra(t)  Ba(t) 


t 

(A-I-C) 

cR*l ,(1) 

(C-A-I) 

cR-l,(l) 

(A-I-C) 

cR=l , (2) 

(A-I-C) 

cR*l ,(2) 

(A-I-C) 

cR=l,(l) 

(A-I-C) 

cR»l.(l) 

(A-I-C) 

cR=l ,(2) 

(C-A-I) 

cr=1,(2) 

0 

0 
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Col: 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

Fig.  6.2 

CR  *  CB  =  ]*  pc  =  ec  =  °-10 
Various  Priorities 

Cases:  Cols.  2  &  6 
Cols.  3  &  7 
Cols.  4  &  8 
Cols.  5  &  9 

(1) :  Strong  Central  Control 

(2) :  Loose  Central  Control 
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PROGRAM  LIFTING  MODEL  I 
(DETER.1INI_.riC  &  STOCHASTIC) 

FILE:  RwAR  FORMA*  A  NAVAL  POSTGRADUATE  SCHOOL 


C 


c 


c 


c 


c 


c 

c 


c 


D 


5 


REAL  ej(  1JOI  ,BA(  IJ')»  ,HC(  1  JOI  »KCC(  1001  »RA(  IOJ) 

HEAL  KC  ,K,h,  THU  A  ,P.U  , HO , RO  ,  RHUUC,  RHjAC 
INTEGER  T  »  T 0 

DATA^TlTLL/Mier^KMlN'.MSTtC  C3','  MODF  L  • ,  •  S  TOCHAST  •  ,  •  I  C  C3 
X  • DEL  '  / 

DATA  1  X/  123436/ »  1P/0/ 

GAVFK  C  J  MHA  r  MO. ILL 
RHUUL=  0 . 

RHOIJA  =  J. 

BUID"  JO. 

BA ( 1 1*0. 

8C 111  =  2  0 . 

RCC (  l)  =  3  0. 

RA( 1 1*120. 

H=  1 . 

T0  =  6  1 


CONT  1IMJF. 

WRITE! 6,  121  I 

READ!  3 .  1  'll )  K  ,KC  .  <  I0AA.RH3AC.  S  ,THJA,  PHI  .  I  TYPE 
IF IK  .LT  .  J. )  STU? 

WRITE! u  « 111  I  K.KC .  UiQJC »  K  H OUA , MHO \A , RHU AC, PH  I ,H, TMUA. S 
WR I T  C ( o  ,  1 1 3 1  ( T 1  I l  J (  I  ,  I T  Y  P  C )  .1=1,3) 

WRITE  I  6, 1)0) 

RO=RCf  t  i  )  *  <  A  I  l  ) 

BO=BA(  1  )  +  HJ  |  1  )  fill.  (  l  ) 

BU  AO=l*  A  (  1)  +  )J(  l) 


DO  10  r=l, TO 
R  1  =  -1  . 

I F I  BUI T I  +  JA ( T )  .Me.  0.)  Rl-RAITI/I  BUITI+dAITI  ) 

R2  =  -l . 

IF  l  BAIT)  .  Mt.  J.  )  K2=K4  I  T  )  /  3AIT) 

R3  =  -  1  . 

IF(RAIl)  .'ML.  .).)  P3*!RA(  T)-RA(in/RA(  II 

ft  4= -  i  # 

I  F  (  T  *.CT.  I  .A. JO.  n  A  (T1. ME.  0.1  K  ,=  (I  A  (  T) -U(  T-l  )  ) /KM  T- 1  I 

WRllTF  (  6  ,  I'll)]  1  l  ,  aCC  (  T)  ,  ,<  A  C  I  )  ,  P  J  1  T  I  ,  :U(  T  )  ,  jC  (  I  )  ,  <  1  ,  R  2  ,  R  3  ,  T  4 

003^*^041  T)  t  ..  J(  T  i  ,  T  00  4  ) 

RCUA* a-J.<  r  ( .:■  i  V  J I  *  i  VJ'.S  l  I  x  ,  I  ,  l  1  vi»r  s 

BUI  T+l  )  =  A  ,  <X  1 1  1  Ul  r  )  -C  I  A'.-  -  JA  ,  ).  I 

dA(  T  +  l  )  *AMi  .1  I  w  A  I  I  )  fC  JA.Of  '  c  J  a  ,  bJ  AO  ) 


RWA00010 
RWA00020 
AWA00030 
KW AO 0040 
MU ' ,RWA0J05J 
RW A0D060 
RwA000  70 
KW AO JO  30 
RWAOOOOO  i 
KWAUOIOO 
RW  A JO 1 10  : 
RW  A  001 20  ; 
RwAOUlAO 
RWAJ01-J 
KWA00150 
RW  AOOlcO 
R.MAJJl  70 
RWAOOiflO 
RW  AOOMO 
Rw A00200 

k*aoo2u  ; 
RWAOJ220  1 
RW  4002  30 
R*AJJ240 
RWA 00250 
RWA00260 
RW AJ027J 
R  W  A  002  d J 
RWA 002*0 
RWA0U303 
RWA 00310 
RW A0032J 
RWA0O330 
Rw AQ0340 
RW 400340  i 
R'W  A003o0  i 
RWA00370  ! 
RWA00330  ! 
*WAO03.?O 
RWAJ0400  I 
R  W  A  004 1 1) 
RW  A0Q42J  ! 
RWA0043J 
RW A00440 
RWA00450 
RWA00460 
RWA00470  i 


cc=ccc  (  rcc  i  r ) » J  •  i ,  :  i  t  j ,  j  j)  . ,  ,  4  , 

SCI  T+l  1  *4  Ux  1 1  .•  v. i  T  )-  it  ♦  3  .  .1  (-L  )  *GA  JSL.Il  X,  IP,  11  YPF  )  ,  >).  ) 

CALL  ALCC  I  Ke  .  I  T  )  ,  JOI  r  I  ,  06  I  n  ,  •  1.I.JC  ,  RH  )«C,M  I  l  ),•'  0,KC,  4C1.AC2) 

RCC ( i  + 1 )  =  \  i.*a i  <  a;citi-im;i  +  c.2I  ► 

X  i.J-  TI  ACi+-.:21  i  GAJSSI  I  X,  IP  ,  i  TYPl  )  ,  J.  ) 

CALL  AL  M  -Ml)  .  J  i<  I  J  ,  lAI  I  j  ,-u;U  ),  RHUi  A.hj,  RHiJAA,  S.C  ,  \L  l  ,  AL  2  ) 


RWA004  3J  i 
**AJJ49J  ! 
RWA00500 
RWA03510 
RW  A0032J 
RWA0053J 
RW  AO 0540 
R  m  A  C05  50 


FILE:  RWAR 


FUHTRAN  A  NAVAL  POSTGRADUATE  SCHOOL 


C 


1 


C 


C 


i 


RAITUI*AMAAl(  HAITI  -  I4LUAL2I  * . 


IS 

l5o 


JFIKA(Tf))«M4|TJ«LT  •  ^•Ul  ^ANO^r!?  }  [U  1  Ucd  |  T  I  !l][  jE  !ol  °An4. 


A  i«nui 

0  CQNT INJt 


TlNJf 

FURMAT  |  /  •  RCCITI  R  4  (  T  )  rtUITJ  B  4 1  T  )  QC(  T  1  »  , 

X  •  RA/dd+3A  R4/RA  «A(  T|-R  A|  d| /UA|  31  RA|  T|  -RA|  W  |/RA|  T- 1 1  '  I 


}?t 


forma 

X  ' 

J2 1  FORMAT 
22  FORMAT 
23  FORMAT 
END 


•1  F^Tfc 
7F'., 


KC  RHQJL  THOUA  RHUA  A  RiiQAC 


PHI' 


.16) 
//JlX, jAd) 


ThJA  S  * . /10F7 .  . 

sK  K,*C  ,RHUA4  ,ftHOAC»S,THJ\,  PHI  .TYPE  (  F6.0M  '  I 


SUBROUTINE  UCCIRCC,  HU, BA  ,RHU'JC,RMUAC,K4,R0,  KC ,  ALC1  ,ALC2| 
REAL  Kw 

ALC1*  RHQJC*RCC*BJ/  Rl7 
A  LC  2*  3. 

JFIRCC  ,LT.  *03001  ,0R,  KC  .L  T.  .  330001)  GO  TO  10 
ALC2*  AH04C^KC*kCC*aA/|KC*RCC  +  R A  ) 

10  RETURN 
END 

SUBROUTINE  ALA(Ra,OU,3A,RC C » KHOJA , RQ  ,KHOAA , KC,  AL 1 . AL2) 
REAL  NC 

AL  1  =  RHCUA  ♦RA*dL'/TO 
AL?=  0. 

IF  IRA  .LI.  .030011  oJ  T'J  10 
AL2*  RHOAA  ml A*3 A/ ( K  CMRCC*  RA ) 

10  RETURN 
END 

FUNCTION  CUAI  IC,  JJ.T  I'JA.M 
RLAL  K 

CUA  =  AM  I  N1  (  T •IJA<’iv<‘au  ,  T HJA«  c)  J  ,  0  J) 

HE  1  UR  1 
END 


FUNCTION  wCi.  lr  CC.OiiJ  ,  1C.B0) 
UCC^  jH1^k  .uVuC/ UJ 
R  t  T  U  < 

END 

F  J  NC  1 1  j-i  .AUSjI  i  »  ,|P,  |  TYPE  ) 
REAL  All  I  10) 

GA J$  S  =  J  . 

IF  I  I  TYPE  ,cj.  1)  RETURN 
IF!  JR  t  JT ,  ))  - J  TO  10 

CALL  L'luM  l  l  A,  1000,1  ,.  ) 

jp=i:oo 

10  GAUSS  =  A  I  i I 


RW AQ0560 
RWA00S7J 
RWA00S30 
KWA0J59J 
RWA00600 
RWA00610 
RWA00620 
RW A  006 30 
RW A00u40 
RWAOJ6SO 
R  WA  00660 
RW AO 06  70 
RW A  006 80 
RW  A  00690 
RWA00730 
RWA0071 ) 
RW A00720 
RW  AO J  730 
R  W  A  00740 
RW  AO J 7 SO 
RWA0076Q 
RWAOU77J 
RW A0078J 
Rw A0J790 
RW A00800 
RW  A  008 10 
Rw A00820 
RWA00300 
RWA00843 
RW  A  QOS  SO 
RWA00860 
R  W  A  008  70 
RW A00680 
RW A  Q0390 
Rw A00400 
RWA03910 
RWA00920 
KW A00930 
RW  A00940 
RW A009S3 
RW A  00960 
RW  A  009  7  0 
RW A  009 80 
RW A 009 90 
RWAUlOOO 
RW AU1010 
RWA01020 
R  W  AO  10 JO 
RW  AO l 040 
RWA01050 
RW AO  1060 
RW  A  0 10  70 
RWA010S3 
RW  A01090 
RW AO l 10Q 


FILE:  RWAR  FORTRAN  A  NAVAL  POSTGRADUATE  SCHOOL 


IP  «  IP 

RETURN 

END 


RWAOll  1(3 
RWA01120 
RWA01130 
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FILE*  WACT 


PH06WM  LISTINQ  MODEL  II  (DETERMINISTIC) 

FUR  TRAN  4  NAVAL  POSTGRADUATE  SCHOOL 


REAL  R ( 3 , 100 1  ,  B  (  3,100)  ,  DR  A  ,Dl)A 

REAL  KK , KB »  RU< 3 ) , BE  (3 )  ,  C HR < 3 ) , CRH  1 3) 

INTEGER  NAMEI3)/ 'A»» '! •»*C •/ 

INTEGER  T  jTO,VKI  3)  ,V0(3) ,  PR  ( 3  ) ,  Pb  (  3  ) 

REAL*8  TITLE!)) 

DATA  TITLE/1 PRIORITY*,*  TARGET T *,' E I)  MODEL'/ 
DATA  IX/12J456/,  IP/J/ 

GAVE«  COMBAT  MUJEL  PRIORITY  TARCETTED 

pji,  u=o. 

R(2,l)=100. 

R(3,  l )  =  20. 

Bii»n*o. 

0(2*11  =  1  bU. 


6(2*11  =  1  bU. 

Bj  3* 11  =  10. 

VR  ( I 1*1=1* 3  ARE  PRIORITIES  FOR  A,1,C  RCSPE.  SAME  FOR  \M!I) 

»m:i 

VR(3  1  =  3 
Vi3 1  l  i  =  1 
VB 12  =2 
VB (33=3 
T0=6l 

CONTINJE 
WRITE! j, 121) 

RE AD (5, 122)  Kft , KB  ,CR.Cfl , (RO(I) *BE(  I  ) *1  =  1,2) 

IF  IKK  .LT.  0.)  STOP 
DO  6  I  =  1  ,  3 
PR!  VK 1  l  )  )  =  I 
PB!  VHj  I )  )  =  I 
CRB! I ) =  -l . 

CBR1 I J»-l. 

CONT INUE 

BAIO=B( l , l ) +3(2,  1) 

PA  1 0=P (1,1) +.3(2,1) 

DR  A  =  - 1  . 

DBA=-I. 

WRITE (6, 111  I  KR.KB.CR.Ca. (ROll  )  ,BE 1  I) ,  I =1,3 
WRITE  !  &,  12  7  )  1  NAME!  VB  I I)  .  I  *1, 3) ,  l  NAME!  /  <U 
WRlTt!3,l23j  (TITLE!  iJ ,1  =  1,3) 

WRITE16, 100) 


KR, KB.C  R 
i NAME  1 VB 
(T ITLt!  I 


CS. (ROll  ) , BE (  l  1,1=1, 3) 

I  I  I  #1  *1, 3)  ,  l  NAME!  /  <!\)  1  ,  i  =  1 ,  31 
,1*1,3) 


DO  1C  T*l,TO 
I  T 1=  T-  l 

WRITE  (6,101)  ITI,(R(I,TJ, 1  =  1*3)  »(•’’(  I ,  T)  ,1  =1  ,  j  )  , 
(CBR!  I )  ,  1  =  1  ,  j)  ,  (CRB i  I  )  ,  I  =  1 , 3  )  , 0 -  A ,  0 ,)A 


CR8(l)=AMlM!  utrjiii.ii.  rt  t  i  ,  i  j  i 

CRB(2)=AMAX1  (0.  ,  AM  INI  (  hi  Pn  (2  ) ,  T  )  ,  k  (  1 ,  l  )-t.  \  1  (  l )  )  ) 

CRB  (3)*A1AXl(0.  ,  AM  1  Nil  B  I  PD 1  3 )  ,  T  )  , R  (  l ,  T  )-wr.  »( 1 )  -  C-:  >  ( 

CBR ( 1 1 =AMl N1 (  K ( PR ( l ) , T ) ,  bll.T)  ) 

CBR (2) =A)AXl(0. ,  AM IMl(R  (PR (2  )  ,T I  , u( 1 , T  )-.  )  <(  l ) )  ) 

CBR  ( 3  )  *  A  MAXI  ( 0.  ,  AM  INI  (  R  ( PR!  3  ) ,  T)  ,  I. !  1 ,  I  )-C  >  •>.  I  1 ) -C  .  s i 


B ( Po I  1 ) 


K  l  l  ,  T ) 


;  i  (  .  )  )  ) 


wacoooio 

WAC0002J 
»ACOOO  3) 
rtA  C  J  00*» J 
.,AC  UOOSJ 
WA  L  0  JO  *»,J 
wAC  00070 
WACO 0 On) 
.,A  C  000°.) 
WACOlllCJ 
wAC  J-Jll » 
OACJOIZ J 
«AC  J  31  30 
W A COO  1 hO 
wACO-Jlil 
•.AC  )  J  1  uJ 
-.AC  )o  17) 
4 A L Jo lb) 
.-.AC  J  II  i) 
-<AcOo2  30 

*  A  C  0  02  ID 
/ .  A  u  J  0  2  2  ) 
0 AC  002 30 
/tACO  j24J 
vif.  002  5  ) 
v  A  C  0  0  2  0 ) 

A  C  J  j  2  T  J 
/.AC  002  OJ 
.\AC0U2DJ 
is  AC0o30D 
.»  A  C  D  0  3 1 0 
.v  A  C  0  J  3  Z  ) 
wA C 003  JO 

/.AC  (3  J  )<•  1 
.JALOJ3bO 
■sAL  3036  3 
. A  w  0  )3  /  J 
*AC.)  33*0 
nAC  0  3330 
s.UJJ“»03 
VAl.  j  J'i  10 
i A  C 00420 
•. AC  J  J  4  3  0 
iAC  J  ■)*,•».) 
sAC  13450 

s  A  t,  J')4  -s.) 

'AC.  00-,  1  j 
•i  A  C  004 J 
..AC  j;;4')o 
..  Al  0  o''J 
•-.Ac  )oS  L J 

•  AC00S2J 
.  ' oOOb  )0 

At.  J  )54J 
.  <i  JMSdO 
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r-  % 


file;  wact 


FORTRAN  A  NAVAL  POSTGRADUATE  SCHOOL 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


19 


ORA»AM 
OBA*AM 
00  XV 
CflRj  I 

CKHii  .  , 

URA*URA-CRb( 
Dd  A«J2A-CBM 
CONTINUE 


NIL  C 

!  3C 
»CbW( 
■CRQI 


CR*KK *»<  (  3f  T  | » CR*RJ  2»  T |  *  R  j  2. T  }  I 
CB*Kb*li|  3  »T  1  »CO*B  J  2*  T I  *  H  1 2 « 1 1  I 


R(  l»T*-l)*AMINX(-<A10i  AN  AX  1(0.,  RU  ,  T  l-C3«  ( VK(  1 ))  ♦  DR4 )  I 
R(2»TUI=*.MINl  (RAIO,  AM  AX l(0»|  K|2'T)-CBR(V4UU  -  CRAI  I 
R|3fTtll*ANlNM  «(i»  U  .  AMAJC  l(UM  R  1 3,TI-CdR  I V  M(.i )  1)1 
B<  l.T  +  X)*ArMU(  iAICf  AlAxKO.,  LI  I  ,  T  l-CAB  U  Jill  I  ♦  OHI  I 
Q  (  2  » T  X I  =  A**  INI  (  0  AI  C  *  AMAX1I0.,  1  (  2  ,  T  I-CRE  I V  0(  2 1  )  -  OJAII 

aijiTni--*<Hi(c;(>,n,  auxiij,,  »(  3,ti-ckjw  ji  u  i  » ) 

XO  CONTINUE 
GJ  TC  5 

< a 1 1 1  *itn  rc  m  xai  ri  aim' , 


100  FORMAT!/* 

X  •  jcm* 


i 


101  FOR1ATCI3foF5.2fJFT.2I 
111  FORMAT  (  /  f '  AR  RJ 

BE (  ROC  BEL*  f 


CR 

/J0F7.31 


CP 


RJA 


ola 


k  :i  i  • 


121  FORMAT  (  *  l  u  ITER  KRf  MAtCR , C  a t RUA ,  dE  A  ,  RO  l  ,  BF  I ,  RJCt  JfcC:  '  I 

122  FOR  -1  AT  (  1  Or  6  .  J ) 


m 


ORMAT  //ilX.i 
FORMAT  '  PKlO 
END 


AB  I 

RI TIES  !  R-BJ 


•  f  3  A2  ,1 2  X  f  *  T-R;  •  ,iA2I 


RAC  00560 
n A  C  J  0  5  70  1 

a AC  0  05  HO 
*AC 00590 
rtACOOrOO 
wA  C  0  Ob  1 0 
mAL00b2J 
,vAC0  0o30  1 

wAC  00640 
* AC00650 
WALOObbO 
*AC  3 Jo  70  ' 

*AC  0 06  00 
mJ  AC  0  06  90 
•<  AL  0 0700  ' 

*AC  )  07  10  . 

.< AC  0  0 720 
■ii  A  l  0  0  7  3  0  ' 
wA C  0  J 7tQ 
.AC  00750 
..AC  00  7  60 
.AC  007  70 
rti*CO  J760 
..Ai:JU79J 
NAC  00800 
m ACOOm 10 
•i  A  L  00  320 
.i  A  C  5  0  3  30  ] 

.1  AC  000  40  ,  ! 
*  \  C  0  'j  d  50  I 
..'ACJJdbO 
..AC  5  JO  70 
.ACOoodO 
w  A  C  0  0  3  90 
•i  AC  009  00 
WAC00910 
ivAC  0 0920 
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